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ABSTRACT 

Some extrasolar giant planets in close orbits—“hot Jupiters”—exhibit larger radii than that of a 
passively cooling planet. The extreme irradiation L eq these hot Jupiters receive from their close in 
stars creates a thick isothermal layer in their envelopes, which slows down their convective cooling, 
allowing them to retain their inflated size for longer. This is yet insufficient to explain the observed 
sizes of the most inflated planets. Some models invoke an additional power source, deposited deep 
in the planet’s envelope. Here we present an analytical model for the cooling of such irradiated, 
and internally heated gas giants. We show that a power source Ldep, deposited at an optical depth 
T<j ep , creates an exterior convective region, between optical depths L eq /L^ ep and r^ep, beyond which 
a thicker isothermal layer exists, which in extreme cases may extend to the center of the planet. 
This convective layer, which occurs only for LdepRiep > L eq , further delays the cooling of the planet. 

Such a planet is equivalent to a planet irradiated with L eq (1 + L^ ep r^ ep /L eq )^, where /? s=s 0.35 is 
an effective power-law index describing the radiative energy density as function of the optical depth 
for a convective planet U oc r fl . Our simple analytical model reproduces the main trends found in 
previous numerical works, and provides an intuitive understanding. We derive scaling laws for the 
cooling rate of the planet, its central temperature, and radius. These scaling laws can be used to 
estimate the effects of tidal or Ohmic dissipation, wind shocks, or any other mechanism involving 
energy deposition, on sizes of hot Jupiters. 

Subject headings: planetary systems — planets and satellites: general 


1. INTRODUCTION 

Hot Jupiters are a class of giant planets that or¬ 
bit their parent stars much closer than Jupiter orbits 
the sun, therefore exhibiting hig h surface temperatures 
(see, e.g. lMavor fc Quelozlll995h . The observed radii of 
many of these exoplanets are greater than Jupiter’s ra¬ 
dius Rj, with the largest planets reaching radii ~ 2 Rj 
(IBaraffe et al.l [20101: Anderson et all 120111: iChan et al.1 


12011b lHartman et al.l 1 20 1 11 ISpiegel fe Burrowsl 1 201 .31) . 

This discovery challenges theoretical evolution mod¬ 
els, which predict that gas giants cool and contract 
to smaller radii, closer to l.ORj, at their inferred 
age, which is often greater than lGyr (IBurrows et al.1 
12000. 120071: [Laughlin et, al.l [2001 . To some extent, the 
larger radii are the result of the intense stellar irradia¬ 
tion, which induces a deep radiative, nearly isothermal, 
layer at the outer edge of the otherwise fully convec - 
tive planets (IGuillot et al.lll995lArras fe Bildst'enll2006l) . 
This radiative layer slows down the evolutionary cool¬ 
ing of the pl anet significantly, in comparison with iso - 
lated planets (IBurrows et al.|[2000i : lChabrier et alIl2004H . 
The slower cooling results in a higher bulk entropy 
at a given age, and theref o re, a larger radius (see , 
e.g. lArr as fe Bildstenl 120061 ISpiegel fe Burrowsl 12012b 
IMarleau fc Cumming|l2014fl . 

Yet, the radii of at least some of the inflated 
hot Jupiters exceed theoretical evolutionary mod¬ 
els, e ven with stellar irradiation taken into ac¬ 
count (IBaraffe et al.ll2003b IBurrows et al.l[2007b iLiu et al.l 
120081) . To settle this discrepancy, a number of 
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The different explanations may be divided into sev¬ 
eral categories. While some works study the effects 
of extra power, deposited at some depth inside the 
pla net, other focus o n enhanced atmospheric opacities 
(IBurrows et al.1 120071) . suppression of convective heat 
loss ( d ouble-diffusive convec t ion, s ee iChabrier fc Baraffel 
120071 iL econte fc Ghabrierl 2012 1, turbulent mixing 
(|Youdin fc Mitchelll 120101 ) ~ or more consistent coupling 
of the day side and night side of the planet, and tak¬ 
ing into account the three-dimensional nature of the 
probl e m (lHansenll2008l lGuillotll2010l ISpiegel fc Burrows 
201C . [20131: Budai et ahl~ 2012b iDobbs-Dixon et al 


2012 : iPerez-Becker fc ShowmarJl2013ilMavne et al.ll2014t 

Rauscher fc Showmanll20l3b 


In this work we focus on the effects of extra 
power sources deep in the planet. Possible heat 
sources inc lude tidal dissipat i on du e to o rbital ec- 
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“thermal tide s” ((Arras fc Soc 
Socratesl I2013D. O h mic h e ating (IBatvgin fc Stevensonl 
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Wu fc Lithwickl 120131 iRogcrs fc Showmanl I2014D . and 
conversion of a portion of the absorbed stellar flux 
into kinetic energy in the a tmosphere, which is then 
dissipated at a greater dep th (IGuillot fc Shownianll2002l 
iShowman fc Guillot 1120021) . While some of these mech¬ 
anisms (e.g. Ohmic heating, kinetic energy dissipation) 
are driven by the stellar irradiation, others (e.g. tidal 
dissipation) are not. If, as some studies suggest, excess 
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planet inflation is cor r elated with the stellar irradiation 
i|Laughlin et al.1120111 : iSchneider et al.11201 111 , then the 
former class of heat deposition mechanisms is favored. 

We do not attempt to determine which power sources 
are more plausible. Alternatively, in the current work we 
study the effects of a general power source on the cooling, 
and therefore the radius, of a hot Jupiter. Specifically, we 
are interested in the dependence of the planet’s cooling 
history on the amount of additional power, and on the 
depth at which this power is deposi ted inside the planet. 
Recently, iSniegel fc Burrowsl (1201,111 have conducted such 
a study numerically (see also lGuillot fc Showmanll2002t 
iBaraffe et al.ll2003l) . In this work we present a simple 
analytical model for the cooling of insolated gas giants 
in the presence of heat deposition somewhere in their en¬ 
velope. Despite our simplifications, this model describes 
the essential processes, and provides an intuitive expla¬ 
nation to the main trends of the scenario. 

The outline of the paper is as follows. In Section [2] we 
present our model for irradiated planets, without addi¬ 
tional heat sources. The effects of additional deposited 
heat are studied in Sections[3]and|4] I n Section[5]we com¬ 
pare o ur results to the calculations offSpicgcl fc Burrow^ 
(j2013T > , and the main conclusions are summarized in Sec¬ 
tion [HI 

2. IRRADIATED PLANETS 

In this section we present our simplified “toy model” 
for irradiated gas giants, without additional heat sources. 
The model is one dimensional, and does not differen¬ 
tiate between_the_day and n ight sides of the planet 
(see ISniegel fe Burrowsl l2013lf . Many aspects of the 
deposition-free model are essentially similar to a previ¬ 
ous analytical work bv I Arras fc Bildsteiil (12006TI , and are 
also implemented by lYoudin fc Mitchelf (|2010t) . 

Isolated gas giants, such as Jupiter, are almost 
fully convective, with convection beg i nning at an op¬ 
tical depth t ~ 1 (iTraftonl I 11964 119671 IHubbard 
119681:[Stevenson fc Salnetedll977l ): see also lGuillot et al. 
(11994) for a reexamination. Insolated planets, on the 
other h and, are expected t o develop a deep radiative en¬ 
velope dGuillot et al.lll996l) . 

2.1. Convective Interior Structure 

The analysis of irradiated planets is simpler when con¬ 
sidering their profiles in the (r, U) plane, with 


[R 

r(r) = / npdr' 
J r 


(i) 


denoting the optical depth at radius r, from the planet 
surface, at radius R , and U = aT A is the radiative en¬ 
ergy density. The density, temperature, and opacity are 
denoted by p, T, and k, respectively, and a is the radia¬ 
tion constant. Our main simplification is to assume that 
this profile can be approximated by a power law in the 
convective interior 


U_ 

TT C 


( 2 ) 


with U c = aT 4 denoting the central radiation energy 
density, determined by the central temperature T c . 

This scaling can be derived in the case of a power- 
law opacity k oc p a T b , as assumed by lArras fc Bildstenl 


(I2006D and lYoudin fc Mitchell (|2010f) . The opacity is 
combined with the adiabatic temperature and pressure 
relation in the convective region (see, e.g., iHubbarcO 
Il977t iSaumon et all [l996h T oc P 1 /( n+1 \ and the re¬ 
sulting polytropic density profile at the edge of the 
planet p « p r ( 1 — r/R) n ( p r denotes th e central den¬ 
sity; see, e.g., iMatzner fc McKeel 1199 9). We neglect 
super-adiabat ic corrections to the convective temper - 
ature profile ( Trafton 19671 iChabrier fc Baraffei 1200 71. 
These power laws yield 


P = 


b+ 1 + n(a + 1) ’ 


(3) 


and t c ~ K c p c R 1 with n c denoting the estimate for the 
central opacity, if the power-law opacity could be contin¬ 
ued to the center. 

The Rosseland opacity (we assume a gray approxi¬ 
mati on) in the range o f inter est is depicted in Figure 
1 of lArras fc Bildsteiil (120061) . The molecular domi¬ 
nated opacity below T 0 « 2500K is fit by a = b = 
0.5. Above T 0 , the opacity rises steeply, with 6 = 7 
(a = 0.5), due to H _ , whose population increases pri- 
maril v as a conse q uence of metal ionization (see, e.g., 
iKippcnhahn fc Weigertl[l994f) . At even higher tempera¬ 
tures, above ~ 10 4 K, the H _ vanishes, due to lack of neu¬ 
tral hydrogen, and the opacity is given by the Kramers 
bound-free opacity law, with a = 1, 6 = —3.5 (see , e.g., 
IKippcnhahn fc Weigertifl994t lPadmanabhanll2000l) . 

As explained in Sections 12.31 and 12.41 the internal lu¬ 
minosity of a planet is determined by the location of the 
radiative-convective boundary. The optical depth (and 
pressure level) of this boundary is given by the opac¬ 
ity power-law which is suitable for the (almost isother¬ 
mal) radiative-zone temperature, and does not depend 
on the deeper opacity structure. For planets without ad¬ 
ditional heating, this temperature is typically b elow Tfi, 
and we therefore get J3 = 0.8, by taking n ~ 2 (IPeeblesI 
11964 IHubbard 119771 : p aurnon et al.11 1996 ). Through¬ 
out the majority of this work we adopt the approxima¬ 
tion of small inflations relative to the zero-temperature 
radius (see Section []}, implying roughly constant cen¬ 
tral densities and pressures. Consequently, a constant 
t c = To ~ 10 11 is suitable for a Jupiter-size planet. In- 
spection of nu meric al evolution models (see Table 1 in 
iBurrows et al.lll997D reveals that the ratio between cen¬ 
tral and photosphere temperatures, which is given by 

setting r ~ 1 in Equation ©, T c /T p h = rf^ 4 , is ~ 100, 
and changes only by a factor of 2 during the cooling of 
a Jupiter mass isolated planet, consistent with our esti¬ 
mate. 

In Section [3] we show that heat deposition raises the 
radiative-convective boundary temperature above To, re¬ 
quiring a more general form, which takes into account the 
strong dependence of the H opacity on the temperature 


II 

to 


(4) 


In this case we get /3 = 0.35 from Equation fl3j), which 
is suitable for temperatures above To. As we show in 
Section l2~3l the convective profile must satisfy /? < 1. 
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2.2. Radiative Envelope Structure 

Proximity to a star dictates an equilibrium temper¬ 
ature T eq on the planet surface, which we take as the 
photosphere, at r ~ 1. Our model is spherical so we as¬ 
sume an even distribution of the a bsorbed heat over t he 
planet surface, for which (see, e.g.. iGuillot et al.lll996f) 


T eq = (1 - A^Tq ' , (5) 

with Tq and R@ denoting the stellar temperature and 
radius, respectively, and D , A are the planet’s orbital 
distance and albedo, respectively. This equilibrium tem¬ 
perature defines an energy density U eq = aT 4 q , and lumi¬ 
nosity L eq = 47ri? 2 crT e 4 q , where a is the Stefan-Boltzmann 
constant. 

The radiation energy is related to the internal luminos¬ 
ity of the planet L; nt through the diffusion approximation 

Lint C dU . . 

4^R2 ~ 3 dF ^ > 

where c is the speed of light. Close to the surface, L; nt is 
constant (in space), since it is the result of cooling of the 
entire planet, and the heat capacity of the outer layers is 
negligible. Therefore, we can integrate Equation ([6]) and 
find the radiative profile 


U = U ( 


3 L h 


eq 


: 4irR 2 


( 7 ) 


where the distinction between the boundary conditions 
U = U eq at r = 0 and r ~ 1 is negli gible, since we 
are ge nerally interested in r > 1. iBodenheimer et ahl 
(120031 use a similar bo undary condition ( with U = 
Upn at r = 2/3), while lArras fc Bildstenl (120061 1 and 
lYoudin fe Mitchelll (I2010D define a temperature TJeep for 
their boundary condition, which is correlated with T eq , 
therefore leading to similar conclusions. More elabo¬ 
rate boundary conditions, which incorporate a frequency- 
dependent opacity, and require the solution of the radia¬ 
tive tra nsfer equation, or the use of a two-stream approxi¬ 
matio n dBarman et al.ll2001l:lHubenv et al.ll2003l : lH anseij 
120081: IQuillotl I2010t IGuillot fc Havell 1201 ill . are beyond 
the scope of this work, and their effect on the temper¬ 
ature deep inside the radiative zone (and therefore on 
the cooling rate) can be modeled by a somewhat differ¬ 
ent effective T eq . We remark that an atmosphere which 
is more transparent in the optical regime (although this 
does not seem to be the case for hot Jupiters; see, e.g., 
iSuiegel fc Burrowsll2013D may induce a greenhouse effect, 
which can be treated with our model as heat deposited 
at r > 1, where the depth for optical light is unity. 

Planar symmetry, and the focus on the outer lay¬ 
ers of the planet throughout this work, are .justified 
both s elf-consistently in Section liOl and bv lGuillot et al.1 
GUI), who find that during the lOGyr evolution of a hot 
Jupiter, the radiative region penetrates only a few per¬ 
cent of the radius. 


2.3. Radiative-Convective Transition 


sets in when 

4 Leg 1 I-P 
3 Lint T < /3 


( 8 ) 


We see from Equation © that at r —>■ 0 the profile is 
radiative. We also see that to ensure convection at r —> 
oo, we must have /3 < 1. 



Fig. 1.— Transition between convective (solid black line) and 
radiative (dashed blue line) profiles. The convective profile follows 
Equation GJ. while the radiative profile is according to Equation 
(0. The combined profile follows the radiative profile for r < r ra d, 
and the convective profile for r > r rac [. Two additional extreme 
cases are plotted for the same values of U c and r c (dotted blue 
lines). A very high value U eq > U c ( 1 — f3) corresponds to a fully 
radiative planet, while a very low value of U e q leads to r ra d < 1, and 
therefore a fully convective planet. Typical values of the optical 
depth and temperature are given. 


The transition between convective and radiative pro- 
files is illustrated schematically in Figure [I] For a planet 
with given interior (C/ c , r c ) and exterior ( U eq ) boundary 
conditions, the convective profile is given by Equation 
0, while the radiative profile is given by Equation 0- 
Since the radiative profile is linear in t, the luminosity 
Li n t and the radiative-convective transition point r ra d are 
calculated by finding the tangent to the adiabatic pro¬ 
file. For r < T ra d, the radiative gradient is smaller than 
the adiabatic, so the profile is radiative. For r > r ra d, 
the radiative gradient is larger than the adiabatic, so the 
profile is convective. Quantitatively, the tangency point 
is found by equating Equations m and 0 and their 
derivatives. The solution for r ra d is 


Trad 

Tc 


p 


1 Ue q 

1 - u c ’ 


which corresponds to a radiation energy density 

G(r ra d) = i _ pU e q , 


( 9 ) 


( 10 ) 


According to the Schwarzschild criterion, when the ra¬ 
diative temperature gradient becomes larger than the 
adiabatic gradient, convective instability develops. Dif¬ 
ferentiating Equations © and 0 shows that convection 


implying that the radiative zone is isothermal to within 
a factor of (1 — /?) -1 / 4 « 1.5. Equation © indi¬ 
cates that the radiative region thickens with increasing 
solar irradiation and with decreasing central tempera- 
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ture (see also IGuillot et all 119961: IBurrows et al.l 120001: 

lArras fe Bildstenll2006l : lYoudin fc Mitchellll2010D . 


2.4. Internal Luminosity and Cooling 

Equations (0) and (0) connect the internal luminosity 
to the central temperature 


Lint ~ ~ 1 (UcV" 

Le q 'Trad Tc \ ^eq / 


( 11 ) 


Note that Equation m breaks down when r ra( j < 1. 
In this case, the planet is considered to be fully convec¬ 
tive, and its luminosity is determined by the radiation en¬ 
ergy density at the photosphere bph- Substituting r = 1 
in Equation © we find 


and luminosity 


U ph = U c 



( 12 ) 


lT = lf L = {7) 0 IT' (13) 

-^eq ^eq \'c/ ^eq 

Comparing Equations © and m shows that the con¬ 
dition r ra d < 1 is equivalent to the condition U p h ^ f7 eq . 
In other words, as long as the central temperature of the 
planet is high enough to ensure that a convective pro¬ 
file will reach r ~ 1 without dropping below T eq , the 
planet remains fully convective. However, as the planet 
cools, the central temperature decreases, and the con¬ 
vective profile would have predicted temperatures below 
T eq . This situation is unphysical, and an approximately 
isothermal radiative zone develops. 

To summarize, using Equations m and m, the in¬ 
ternal luminosity of the planet, as a function of its central 
radiation energy density is given by 



■Hn. < T P 

? 

(14) 

LLl. > 

Req ^ T C 


~ -^eq /-^int n 

^ 1. Since 


where the transition is when r rad 
/3 < 1, Equation m indicates that the proximity to 
a star, which is manifested for U c /U eq < , dimin¬ 

ishes the internal lu minosity and s l ows d o wn the planet’s 
cooling (see also IGuillot et al.1 119961: IBurrows et al.1 


20001: lArras fc Bile sten 120061 lYoudin fc Mitchelll 120101 
Sniegel fc Burrowsll201,lH . 

As a planet cools, it makes the transition from an iso¬ 
lated planet, which is not affected by stellar irradiation, 
and matches the high U c limit of Equation (fT4l) . to an 
insolated planet, which matches the low U c limit. The 
age, at which this transition takes place, may be esti¬ 
mated by writing an evolution equation for the central 
temperature of an isolated planet 


M dT c 
m p dt 1 


(15) 


where ks is Boltzmann’s constant, m p is the proton 
mass, M is the mass of the planet, and t denotes 
time. Equation m is valid for the cooling degen¬ 
erate phase of the planet’s evolution, which follows 


a rapid non-degenerate cont raction phase (see IGuillot! 
I2005t lArras fc Bildstenll2006D . Substituting the luminos¬ 
ity from Equation (O, we get 

4tt R 2 aT* 


Tc 


M_dTc 
B m p dt 


(16) 


Assuming that the planet cools down from high tem¬ 
peratures, Equation (fl6l) predicts the cooling time to a 

temperature T c = T eq Tc ^ 4 , which marks the transition 
to the insolated regime (at this stage the photosphere 
reaches a temperature T p h ~ T eq ): 


t 


M k B 0/4 
m p 127Ti? 2 crT e 3 q c 


7Myr 



(17) 


were the numerical value is estimated for Mj (Jupiter 
mass) and Rj. T his estimate r oughly fits Jupiter, which 
is a few Gyr old (lGuillotll2005l ). has an equilibrium tem¬ 
perature T eq ~ 11 O K, and is close to the transition, 
since T p h « 120K (jHanel et all 1 198 il l Strongly irradi¬ 
ated planets orbit their parent star closer than Jupiter 
by a factor of 100, thus having, using Equation ©, 
T eq ~ 10 3 K, and therefore reach the insolated regime 
and develop a r adiative zone after less than 10 7 years 
(consistent with IGuillot et all 1 19961) . 

In the following, insolated stage of planetary evolution, 
Equation (1151) is combined with the low U c limit of Equa¬ 
tion (fl4l) to produce a cooling equation 


Tc(t) 


M T eq (1 “ /3)//3 1 
c B m p AnR 2 cr t 

1/4 


/3/(4-/3) 




1/4 


(18) 


with the numerical valu es estimated f or Mj , Rj. This 
result is consistent with IGuillot et all (I1996H . At an age 
of 8Gyr, the radiative region encompasses a few percent 
i n radius, 1 — r/R~ (p/p c Y^ n ~ T eq /T c , consistent with 
IGuillot et al.1 (119961) , and self-consistently justifying our 
focus on the outer layers of the planet. 

A third stage in a cooling planet’s life, equilibrium, 
is reached when U c = U eq . In this case, according to 
Equation ©, the planet becomes fully radiative, and the 
luminosity vanishes Lj n t = 0 (since the temperature gra¬ 
dient vanishes). This condition determines the planet’s 
equilibrium central temperature, and therefore its equi¬ 
librium radius. More precisely, the planet becomes fully 
radiative when U c drops below U eq /(1 — (3). It is easy to 
see, from Equation (fl8l) . that without additional heat de¬ 
position, planets never reach equilibrium (the timescales 
are too long). 

The three stages of planetary evolution, described in 
this section, are displayed in Figured] for planets with 
the same central temperature and optical depth, but dif¬ 
ferent equilibrium temperatures. 


3. DEPOSITION OF ADDITIONAL POWER 


In this section we discuss how the internal luminosity 
is affected by an additional heat source, which deposits 
power at some depth inside the planet. Our goal is to 
find L in t(T eq ,T c ,L dep ,r d ep), with L dep denoting the de¬ 
posited luminosity, and r dep denoting the optical depth 
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where the power is deposited. In Section 12.41 and specif¬ 
ically in Equation (1141) . we found L; nt (Xdep = 0), which 
we abbreviate hereafter as L; nt . We focus on plan¬ 
ets which have developed a substantial radiative zone, 
Trad > 1, so L? nt is determined by the low U c limit of 
Equation (1141) . The motivation is that strongly irradi¬ 
ated planets develop such a radiative zone in less than 10 7 
years , as explained in Section l2~4l isee also lGuillot et all 

SMI). 

Without heat sources, the radiative zone is described 
by Equation ©. The addition of a deposited power Ldep 
at Tdep introduces a jump in the profile slope. Explicitly, 
the radiative profile is altered: 


dU _ 3 

dr AttR 2 c 


Ltot — Lint T Tdep T <C Tdep 


L 


int 


T > Tdep 


, ( 19 ) 


with Ltot denoting the total luminosity for t < Tdep- 
The radiative profile, with deposited heat included, is 
displayed in Figure [2] As seen in Figure [21 the deposited 
heat effectively changes U eq to some higher U® ® . T he 
internal luminosity Li nt is found, as in Section 12.31 by 
calculating the tangent from the point [Td ep , U(Td ep )] to 
the convective profile, which is equivalent to calculating 
the tangent from the point [0,17®®]. Therefore, the re¬ 
sults of Section [2731 are reproduced, but with U/r® instead 

of U eq . Since, from Equation 0ED- L int OC Ueq 1 ^'' P , the 
internal luminosity decreases. This result is also evident 
graphically in Figure [2] the slope of the radiative tan¬ 
gent decreases, as well as the deeper penetration of the 
radiative zone into the planet interior. 


3.1. Secondary Convective Region 

In Figure [2] we assumed that the outer layers remain 
radiative. However, if the total power, L to t, which has to 
be evicted from the planet, is high enough, then it cannot 
be delivered radiatively, and a secondary convective in¬ 
stability develops. According to Equations 0 and 03), 
with L t ot replacing Li nt as the evicted power, convection 
appears at an optical depth 


4 /3 Leg 
3 1-/3 L tot 


( 20 ) 


where the radiation energy density reaches 17 eq /(l — /3). 
If we focus on intense deposition Ldep 7S> Li n t, then we 
may approximate L tot ~ Ld ep , and using Equation (I2U1) . 
Tb ^ L e q/Ldep- 

The temperature profile, with deposited heat added, 
can be divided into two cases, depending on Ld ep and 
Tdep- If Tdep < T b , or equivalently L d epTde P /L eq < 1, The 
depiction in Figure [2] is accurate. As seen in Figure [2j 
in this case C/™ is restricted by a maximal value f7|® < 
f7 eq /(l— /3), limiting the reduction in internal luminosity, 
using Equation ED. to an insignificant effect L m t /LP nt > 
(1 — „ l. 

In the case where Td ep > Tf,, a secondary convective 
region emerges between Tb and Td ep . This region is dis¬ 
played in Figure [3] The modified luminosity can be cal¬ 
culated by finding t/|®. Using Equation (l20l) . approx¬ 
imating Ldep O’ Li n t, and noting that in this regime 



Fig. 2.— Modification of the radiative profile due to the depo¬ 
sition of heat at optical depth Tdep- The convective profile (solid 
black line), and the radiative profile without a heat source (dashed 
blue line) are similar to Figure [[] The modifie d rad iative profile 
(dot-dashed red line) is described by Equation (1191 . The imagi¬ 
nary continuation of the profile with a slope 3Li nt /(47ri? 2 c) down 
to r = 0 (dotted red line) defines . As in Figure [T] the com¬ 
bined profile follows the radiative profile up to the tangency point, 
and then it follows the convective profile. 



Fig. 3.— Emergence of a secondary convective region from 
to Tdep (solid green line). The convective (solid black line) and 
radiative (dashed blue line) profiles without a heat source are sim¬ 
ilar to Figures [T1 and [2l The radiative profile with heat deposited 
at Tdep (upper dot-dashed red line) is similar to Figure [2] and is 
inadequate in this case. Due to the secondary convective region, 
the radiative profile is altered, and it follows the lower dot-dashed 
red line, up to the tangency point with the main convective profile. 


-^ 7 depTdep/-^- / eq ^ lj WG find 

^eq 1 U(T dep) / L dep T dep \ P 

U eq ~ 1 - /3 U(T b ) T L eq j ’ ^ j 

where U(Td ep ) is adiabatically related to U{rb). Com¬ 
bining this result with Equation (fTT1) . we get Lj n t/L? nt ~ 
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(LdepTdep/ieq) ^ which we rewrite as 


Tint 




TdepTde 


T. 


eq 




( 22 ) 


in order to interpolate with the L dep T dep /L eq 1 regime. 

Equation C^l) shows that the figure of merit of the heat 
source, in terms of its influence on the internal luminos¬ 
ity, is LdepTdep/ieq, which must satisfy L dep T d e P /L eq > 1 
for a significant effect. For example, if the heat is de¬ 
posited at Tdep ~ 1, then T<j ep > T eq is needed. This 
result is intuitive, since the equilibrium luminosity it¬ 
self is deposited at r ~ 1. If, on the other hand, the 
heat is deposited at the top of the convective region, 
^dep ~ T rad , then, using Equation (HI]), T dep > L? t 
is sufficient. Although not evident from Equation (1221) . 
Lden > is required fo r a significant effect (see also 
IGuillot fc Showman! 120021) . even for deposition in the 
convective region Td ep > T ra d- Otherwise, the deposited 
power simply replaces part of the internal luminosity as 
the convectively delivered heat which maintains the ra¬ 
diative profile at the envelope, and the temperature pro¬ 
file remains unchanged (with respect to Tdep = 0). In 
this case, the luminosity decreases by Tdep- 

The secondary convective region, described in 
this section, was als o no ted in numerica l works 
JGuilloH^Showman I2002t IBatvein et all 120111 : 
IWu fe Lithwickll2013 ). for high dissipated power. 


3.2. Cooling and Equilibrium 

Substantial additional heating raises the effective equi¬ 
librium temperature, and therefore the inner radiative- 
convective boundary temperature above To, as seen in 
Equation (T2T1) . Under the influence of such a heat source, 
using Equation ([221) . the cooling Equation (fTHl) in the in¬ 
stated regime changes to 


T c {t) oc 



TdepTdi 


ep 



TdepTdi 


ep 


Jeq 


/3(l-/3)/(4-/3-/3b) 

0.19 


(23) 


slowing down the planet’s cooling, with the power calcu¬ 
lated using the H _ opacity values /? = 0.35 and 6=7, 
following the discussion in Section [2] 

Equilibrium is reached when U c « U^ q (see Section 
[2|. According to Equations (0 and (12T1) . the required 
heating intensity for this scenario is 

TdepTdep f Cc 

Teq \ U e q 



with r rac j denoting, as before, the original radiative- 
convective transition, without additional heating. We 
may rewrite Equation < 0 , using Equation CEB , as 
TdepTdep ~ TP nt r c , showing that, as intuitively expected, 
a deposition of order TP nt forces equilibrium, when the 
deposition is at the center. For a given heat source, 
Equation (1241) defines the equilibrium central temper¬ 
ature T c = T®® = (Ul^/a) 1 / 4 , for which evolutionary 
cooling not only slows down, but stops entirely 


IE = 1*1 „ (l , ^epTdep 

T eq T eq \ T eq J 


( 25 ) 


where we have again interpolated with the weak heating 
regime. 

If, due to very intensive (or deep) heating, T®® > 
10 4 K, the opacity behaves according to the Kramers 
law, with b = —3.5 and (3 > 1, prohibiting convec¬ 
tion (see Section [B • In this case, the planet becomes 
fully radiative from Td ep inward, even before equilib¬ 
rium is reached. The internal luminosity in this case 
is given by Equation 0, T int ~ (4ttR 2 c/3)(U c /t c ). 
Combining this result with Equation (|15l) . we find that 

the cooling time is t oc T c ^ = T c -6 ' 5 , with hotter 
(more inflated) planets reaching equilibrium first (plan¬ 
ets with radii close to 2 Rj cool even faster, due to 
the dependence of the opacity on the density). This 
regime corresponds to the equilibrium-inflation paradigm 
of hot Jupiters, which typically involves very deep 
heating ( IGuillot fe Showmanll2002t iChabrier et al.l12 001: 

Burro ws et al 1 120071 iLiu et al.ll2008t iSniegcl fc Burrowsl 

2013). 

4. EFFECT OF HEATING ON PLANET RADIUS 

In Sections[2]and[3]we found the effects of stellar irradi¬ 
ation and additional heat deposition on the planet lumi¬ 
nosity. The diminished luminosity results in slower cool¬ 
ing, and therefore a higher bulk temperature and larger 
radius at a specific age, in addition to a larger equilib¬ 
rium radius. For a rough estimate of this inflation, we 
introduce a simple model which relates the plant’s radius 
R to its central temperature T c . 

We start by writing a simple equation of state (EOS) 


P = Kp 5/3 



(26) 


where K ~ h 2 /(m e m^ 3 ) is the electron degeneracy 
term, and ( po/m p ) 1//3 ~ m e e 2 /h 2 ~ a^ 1 is the elec¬ 
trostatic correction ( h denotes the Planck constant, m e 
and e denote the electron mass and charge, respect ively, 
and a 0 is the Bohr radius; see, e.g.. !Padmanabhanll2001l 
where we assume a mixture of hydrogen and helium, 
so the atomic weight and charge are approximately 1). 
The rightmost term is the ideal gas pressure of the non¬ 
degenerate ions. We are interested in the regime where 
the degeneracy parameter 6 = (pksT / m p ) / Kp 5 ? 3 ~ 
HU 1 , which is adequate for HD 209458b (see iGuillotl 
12005 ). In this regime, the electrons are almost completely 
degenerate, and the ion thermal term is a small correc¬ 
tion. Since the degeneracy pressure is inversely related to 
the particle mass, the ions are non-degenerate. The ther¬ 
mal contributio n of the degenerate ele ctrons is quadratic 
in 9 (sec, e.g.. I Ghandrasekhail 119391 ) and is therefore 
neglig ible. This EOS is also used by lArras fc Bildstenl 
(120061) . 

We proceed by replacing the profiles p(r),T(r), P(r), 
with the characteristic values p c ,T c ,P c . This simplifi¬ 
cation is justified by noting that the entire planet (ex¬ 
cept for the radiative isothermal envelope) can be de¬ 
scribed approximately with a single almost degenerate 
polytrope. The distinction between a degenerate core 
and a non-degenerate envelope is insignificant, since the 
te mperature pr ofile is almost parallel to lines of constant 
9 (lGuillotll2005D . A similar argument renders the distinc¬ 
tion between a core without electrostatic corrections, and 
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an envelope with significant corrections, also insignificant 
in comparison with the bulk electrostatic correction. For 
a given central temperature T c , the contribution of the 
isothermal (T « T eq ) radiative envelope to the radius 
is negligible, because it scales with the temperature, and 
T eq <C T c . We emphasize that the quasi-degenerate treat¬ 
ment in this section only holds for mild inflations, and 
should be taken only as an order of magnitude estimate 
for extremely inflated planets, with radii ~ 2 Rj. In these 
cases, numerical radius-temperature relations are more 
adequate, as detailed below. 

In equilibrium, the pressure of Equation (l26l) balances 
the gravitational pressure, of order GM 2 /Ft 4 (G is the 
gravitational constant), 


Kpl 1 3 


1-1^ 


1/3' 


+ —k B T c ~ GM 2 / 3 p 4 J 3 . (27) 
m p 


For small variations of the temperature, differentiating 
Equation (1271) . while holding the planet mass constant, 
yields a linear relation 


dR 

~R 


M 


1 dpc 

3 p c 


—k B dT c 

m __ ■ LJ 


M 


KpT - b t c 


(28) 


In the degenerate limit ( p c k B T c /m p <C Kp 3 ^ 3 ), the 
radius increases with temperature, while in the non¬ 
degenerate limit ( p c k B T c /m p Kp 3 J 3 ), the planet 

shrinks as the temperature increases. As mentioned 
above, the non-degenerate limit corresponds to the initial 
contraction phase, accompanied by an increase in tem¬ 
perature, while the degenerate limit describes the later 
cooling of the planet, which is accompanied by a decrease 
in radius. The inflated hot Jupiters are on the late-time 
cooling branch of this scenario (the tran sition is at a ra¬ 
dius around 4 Rj, see lGuillot ct al.lll996f) . 

We make a linear approximation for low temperatures, 
using Equations (1771) and (1751) 


A R. 


R 2 


GMm 


-k B T c 


1 + ( Mj/M ) 213 


O.L R 


\10 4 K J 


k B T c 

m p g 


(29) 


with A R = R — i?o denoting the inflation relative to 
the zero-temperature radius Rq, g is the surface gravity, 
and noting that Jupiter’s mass is close the well-known 
inversion of the zero-temperature mass-radius relation 

Mj ~ (K/G) 3 I 2 p^ 2 dPadmanabhanl [200111 . Since we 
focus on hot Jupiters with M ss Mj , the electrostatic 
correction term Mj/M does not influence our order of 
magnitude estimate. This relation fits well numerical 
radiu s-temperature curves (for example, iBurrows et al.l 
mu, Table 1), up to a radius of ~ 1.3i?j, where 
the linear approximation breaks down (sin c e it as sumes 
A R <C i?o ~ Rj). See lArras fc Bildstenl (|200fiD for a 
similar linear relation between the radius change and the 
central temperature. 

Although Equation (1771) was derived for a constant 
planet mass, it is also valid (up to numerical factors) 
when the planet surface gravity g , instead of the planet 
mass is kept constant. In this case, the gravitational 


pressure in Equation (1271) may be written as g 2 /G, and 
the radius increase is found by differentiating dR/R\ g = 
—dp c /p c \ g . This s cenario is interesting for comparison 
with the results of lSoiegcl fc Burrowsl (12013 1. who stud¬ 
ied models with constant surface gravity. 


5. COMPARISON TO NUMERICAL CALCULATIONS 

iGuillot fc Showmaiil (120021) and iBaraffc et al.l (120031 ) 
find that less heat dissipation is needed to inflate hot 
Jupiters if the heat is deposited in deeper m ass fractions 
of the planet (see also lWu fc Lithwickll2013D . In this sec- 
tion we compare our estim ates to the systematic work of 
iSmegel fc Burrowsl (120131) . henceforth SB2013, who stud¬ 
ied the effects of different amounts of additional heat, 
deposited at different pressure levels. 

In their Section 3.3, specifically Figure 4, SB2013 de¬ 
posit different fractions of the incident flux Fg ep /F eq = 
Adep /L eq at different pressures inside a planet similar to 
HD 209458b. For each value of additional power and 
deposition pressure, SB2013 present the inflated planet 
radius. We emphasize that, in their Section 3.3, SB2013 
concentrate on planets which have not yet reached equi¬ 
librium, but are still contracting. Moreover, they com¬ 
pare models with the same effective temperature, defined 
as T e g = (Tint/oTTri? 2 ) 1 / 4 , and the same surface gravity, 
resulting in models with different ages and masses. 

By combining Equations (1171) and (1771) . we relate the 
central temperature to the deposited power, 


Tc_ 

T eq 


( T eq y / TdepUdep 

Tq / T e q V T eq 


1-/31 


/3/(4—0 b) 


(30) 

The inflated planet radius can be related to the deposited 
power by combining Equation (1301) with Equation (1291) 


A R = 


k B T e q 

m p g 


TO 


TeS \ 4 ( T eq ' 61 N{i ^ 
T 


eq / \ Tq 


X 1 + 


TdepUde 


/3(l-/3)/(4-/36) 


(31) 


= Ai? 0 (l + LdepTdep 

T e q 


/3(l-/3)/(4-/36) 


with Ai? 0 = 0.27 Rj denoting the expansion without ex¬ 
tra heating. This numeric value is estimated by substi¬ 
tuting the parameters g = 10 4 cm s -2 and T e g = 180K, 
which arc kept constant in SB2 013, the observed tem - 
perature T eq = 1300JF T e g (|Crossfield et al.ir2012lh 
and taking the adequate /? = 0.35 and b = 7. The 
model of SB2013 for HD209458b is indeed inflated by 
AR 0 = 1.25i?j—0.95i?j = 0.30i?j, due to the irradiation 
of t he parent star, even without extra power s ources (see 
also iBaraffe et al.l 120031 : IBurrows et al.l 120071) . in agree¬ 
ment with our estimate. The zero -temper a ture r adius 
Ro = 0.95/7/ is e stimated using iGuillotl (|2005l) and 
IBaraffe et al.1 (120101 ) , for the given surface gravity. An al¬ 
ternative estimate of the inflation may be given by relat¬ 
ing the central temperature, found in Equation (1301) . to 
A lt, using the numerica lly calculated R{T C ) curve given 
in IBurrows et all (jl997l) . This analysis is more accurate 
at large expansions. 




































































In Figure 0] we compare the results of SB2013 (their 
Figure 4) to the equation 


A-R d ep — Ai? 0 


I +77 


TdepAie 


- 1 


(32) 


with ARdep = A R — ARq denoting the expansion rela¬ 
tive to the radius without extra heat deposition (1.25I?j 
in our case), A R 0 = 0.30 Rj, 77 an order unity fitting 
parameter, and <5 = /3(1 — /3)/(4 — /3b) ~ 0.15 from Equa¬ 
tion (f3Tj> . The fraction of incident flux Ld ep /L eq is given 
by SB2013, and the translation of deposition pressure to 
optical depth is done approximately b y P/g «t/k, with 
k taken from I Arras fe Bildstenl (120061) . and uncertainties 
going into the fitting parameter 77 . The parameters in 
Figure O 77 = S = 0.2, are not, strictly speaking, a best 
fit, but only an illustration. A second fit, with 77 = 0.25 
and the predicted S = 0.15, is also provided. This fit, 
however, incorp orates the more acc urate nonlinear ra¬ 
dius correction dBurrows et al.l[T997h to Equation (l32l) . 

While our model is not exact, as seen in Figure 0] it 
does manage to describe the main trends of atmospheric 
heating, and, given the general uncertainties and ap¬ 
proximations of the model, it fits the numerical calcu¬ 
lations reasonably well over several orders of magnitude 
in pressure level and flux. When the realistic radius- 
temperature relation is used, the analytically predicted 
power 5 = 0.15 fits the results adequately. With the ap¬ 
proximated linear radius-temperature relation, however, 
we fit S = 0.2 to the numerical results. 



Fraction of Incident Flux 


Fig. 4.— Comparison of our prediction for planet inflation (rel¬ 
ative to 1.257?./), due to heating (expressed as a fraction of the 
incident flux) at different pressure levels, with lSpicgel fc Burrow j 
( 121 ) 1 -il l, for constant surface gravity g = 10 3 cm s~ 2 , effective tem¬ 
perature T e fj = 180K, and incident irradiation of HD 209 458b . Our 
analytical model (dashed blue lines) is given by Equation (1321) . with 
fitted parameters rj = S = 0.2. A second model (dot-dashed red 
lines), which takes into account the nonlinear radius-temperature 
relation, is also provided, with a fitted 77 = 0.25, and the analyt¬ 
ically derived <5 = 0.15. The results of SB2013 (solid black lines, 
marked with circles) are taken from their Figure 4. 


In the “weak heating” regime TdepTdep/Teq < 1, it is 


possible to linearize Equation (13^1) 

AA dep « S V AR 0 = (AR 0 , (33) 

-^eq T/ e q 

which has the advantage of a single fitting parameter £ = 
Sr/ ~ 0.04, allowing for a more intuitive interpretation of 
the results. 


6 . CONCLUSIONS 


In this work we presented a simplified analytical model 
for irradiated giant gas planets, which includes an addi¬ 
tional heat source, deposited at an arbitrary depth inside 
the planet’s atmosphere. 

For an irradiated planet with no extra heat sources, we 
found useful approximate scaling relations for the inter¬ 
nal luminosity L lnt of the planet, and for the radiative- 
convective boundary, marked by an optical depth r ra( j. 
Explicitly, the radiative-convective boundary and the in¬ 
ternal luminosity are given by 



(34) 


with L eq , U eq oc T e 4 q denoting the equilibrium luminos¬ 
ity and radiation energy density imposed by the stellar 
irradiation, respectively (T eq is the equilibrium temper¬ 
ature). U c is the central radiation energy density of the 
planet, and r c represents its optical depth (with opac¬ 
ities extrapolated from the outer layers of the planet 
inward). The power /3 ~ 0.35 is related to the adia¬ 
batic index of the planet and to the opacity near the 
planet’s edge. In addition, we identified the condition 
Trad ~ 1 , as the boundary separating insolated (irradi¬ 
ated) planets (r ra d > 1 ) from isolated ones (r rad < 1 ). 


well-known results (IGuil 

ot et a 

. 1996|: Burrows et al. 

2000 

: Arras & Bildsten 

20061 

Fortncv & Nettclmann 

201C 

;lYoudin & Mitchell 201(1 ISoiegel & Burrows 20131). 


We addressed the effects of possible additional heat 
sources, parametrized by additional luminosity L d ep, de¬ 
posited at an optical depth Td ep . Such heat sources lower 
the internal luminosity of the planet, slow its evolution¬ 
ary cooling, and are therefore candidates f or sol ving the 
puzzl e of over-inflated hot J upiter s (see Baraffe et al 


201C; iFortnev fe Nettelmannl 120101 : iSoiegel fe Burrows 
20ll iBaraffe et all 12014 and references within). Ac¬ 


cording to our model, the figure of merit in this context 
is LdepTdep/Aeq, with idepT de p/Teq > 1 necessary for a 
significant effect. Concretely, an outer convective layer 
is formed in this case, between optical depths T e q/T d ep 
and Tdep, which is equivalent to an enhanced stellar irra¬ 
diation 


t eff 
-^eq 


'yeff ' 
1 eq 

Tea 


1 + 


TdepTde 


(35) 


The internal luminosity is therefore reduced from a value 
of L? nt without extra heating to, roughly, 


Mnt 


Y , -bdepTdep ) 
Teq / 


-(1-/3) 


(36) 


with a further threshold of L® nt required for a significant 
effect. 
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Since the inflation of a planet, relative to its 
zero-temperature radius, is proportional to its cen¬ 
tral temperature, these relations provide an in¬ 
tuitive explanation to the dependence of in f lation 
on the deposition depth (I Guillot fe Showmanl 120021 : 
IBaraffe et al.1 120031 IWu &; Lithwickl 120131! . and to the 


heating (IGuillot & Showmanl 2002 

Burrows et al.l 20071 

Liu et al.ll2008l ISuiegel & Burrows 

1201311. 


Quantitatively, we compared our mode l to th e recent 
parameter survey of iSpiegel fo Burrowsl ( 201 3!). The 
model, though not exact, fits the numerical results ade¬ 
quately, considering its simplifying assumptions and ap¬ 
proximations. 


The scaling laws presented in this work, combined 
with potential observational correlation s (see, e.g., 
iLaughlin et al.ll20TTI : iSchneider et al.l[2011ll . may be used 
to rule out or to favor specific energy sources, or even the 
entire additional heat source paradigm. 
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